Cold nuclear matter physics at forward rapidities 
from tf+Au collisions in PHENIX 
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00 | Abstract. We present measurements by the PHENIX experiment at RHIC of di- 

0^ ' hadron pair production in d+Au collisions where the particles in the pair are varied 

across a wide range of pseudorapidity, out to T) = 3.8. With di-hadrons, varying 
£><J I the pt and rapidity of the particles in the di-hadron pair allows studying any effects 

as a function of partonic x in the nucleus. These di-hadron measurements might 
probe down to parton momentum fractions x ~ 10 -3 in the gold nucleus, where the 
' interesting possibility of observing gluon saturation effects at RHIC is the greatest. Our 

measurements show that the correlated yield of back-to-back pairs in d+Au collisions 
is suppressed by up to an order of magnitude relative to p+p collisions, and increases 
with greater nuclear path thickness and with a selection for lower x in the Au nucleus. 



PACS numbers: 21.65. Qr,25.75.-q,25.75.Bh 



ON 

1. Introduction 



Deuteron-gold collisions at RHIC provide a means to explore nuclear effects on the 
initial-state parton densities in the nucleus, which is vitally important for understanding 
the baseline production in heavy-ion collisions. RHIC experiments have shown that 
single inclusive hadron yields in the forward (deuteron) rapidity direction for y/s NN = 
200 GeV d+Au collisions are suppressed relative to p+p collisions [TJ [2j E] . The 
mechanism for the suppression has not been firmly established. Many effects have 
been proposed for this suppression, such as gluon saturation [U [5] , initial state energy 
loss [6j [7], parton recombination [8], multi-parton interactions [9], and leading and 
higher-twist shadowing [TOj [11] . 

One set of measurements that might help to distinguish between the competing 
models is forward azimuthally correlated di-hadron correlation functions, which directly 
probe di-jet production through their 2—^2 back-to-back peak at A0 = 7r. This technique 
has been used extensively at RHIC and is described in detail elsewhere [121 [El E]- 
The di-hadron results presented here were obtained from p+p and d+Au runs in 2008 
with the PHENIX detector and include a new electromagnetic calorimeter, the Muon 
Piston Calorimeter (MPC), with an acceptance of 3.1 < 77 < 3.8 in pseudorapidity and 
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< 4> < 2tt. Di-hadron measurements can probe more precise ranges of parton x in a 
gold nucleus than do single hadron probes (e.g., R&a). At forward rapidities, a single 
hadron probe will cover a very broad range of x, 10~ 3 < xau < 0.5, thus mixing together 
the shadowing, anti-shadowing, and even EMC effects [10]. Azimuthally correlated di- 
hadron measurements also enhance the di-jet fraction in the event selection, since one 
selects only the back-to-back hadrons. 
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Figure 1. Background-subtracted mid- forward rapidity 7r°-7r° (top) and forward- 
forward cluster- 7T° (bottom) per-trigger correlation functions, for p+p (open points), 
d+Au peripheral (60-88%, triangles), and d+Au central collisions (0-20%, solid points) 
at y/sNN = 200 GeV. The trigger px ranges from 1.1 to 7 GeV/c and the associated 
7r°'s have pr = 0.5 — 0.75 GeV/c. Systematic errors of up to 30% on the near side 
{\A(p\ < 0.5) have not been shown. The subtracted pedestals, bo, are shown for each 
case. 



By performing several correlation measurements with particles at different pr and 
rapidities, one can systematically scan different x ranges with an observable that is 
enhanced for the leading-order perturbative QCD component. Probing the x dependence 
of the effect is an important test since most models predict that any effects should be 
stronger at smaller x. Particles at higher pseudorapidities are produced from smaller x, 
so measuring hadrons from more forward rapidities should probe smaller x. 

2. PHENIX MPC d+Au di-Hadron Correlations 

For this analysis, back-to-back 7r°-7r or hadron- 7r° pairs are measured with one particle 
at mid-rapidity, and the other at forward rapidity. Back-to-back cluster- 7r° pairs are also 
measured where both are in the forward rapidity region. The clusters are reconstructed 
from the energy deposit of photons in the MPC, and are estimated to be at least 80% 



Forward di-hadron measurements in d+Au collisions with PHENIX 



3 



dominated by 7r°'s, with the remainder coming from single photons from decays of 77's 
and from direct photons. Further details of the analysis are available in |14j . 

As shown in figure [TJ the away-side peak for d+Au central collisions appears 
significantly suppressed compared to p+p collisions and peripheral d+Au collisions. 
This effect is large for the mid-forward di-hadron correlations and becomes even larger 
for the forward-forward correlations. Within large errors, the Gaussian widths of the 
away-side correlation peak for the mid-forward di-hadron correlations remain the same 
between p+p and central d+Au. For the forward-forward case, uncertainties in the 
pedestal level from the underlying event and the strong suppression of the away-side 
peak make extracting the width unreliable. For this case, the away side peak width in 
central d+Au collisions is allowed to vary up to twice as much as in p+p when accounting 
for this systematic uncertainty. 




Figure 2. JdA versus x A ™ a for peripheral (60-88%, in red) and central (0-20%, in 
blue) d+Au collisions at = 200 GeV, and compared to the EPS09 LO Rf u (xAu) 

curves at a scale Q 2 = 4 GeV 2 [15] . 



The observed suppression is quantified by studying the relative yield, JdA [IS], of 
correlated back-to-back hadron pairs in d+Au collision compared to p+p collisions scaled 
with the average number of binary nucleon collisions (N co u), 

1 pair I 

j 1 a dA / a dA m 

dA / AT \ P air I ' ^ ' 

\A coll) CTpp /Opp 

This is simply the analog of the usual nuclear modification factor i?dA but for hadron 
pairs. The a are the p+p or d+Au inelastic cross-sections, while a pair is the cross-section 
for di-hadron pair production. The JdA is calculated using the correlated away side 
peak after subtracting the pedestal 60 • ^dA decreases with increasing number of binary 
collisions, (N co u), or equivalently with increasing nuclear thickness. The suppression 
also increases with decreasing particle pt and is significantly larger for forward-forward 
hadron pairs than for mid-forward pairs. The observed suppression of JdA versus nuclear 
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thickness, p? and r] points to large cold nuclear matter effects arising at low parton 
momentum fractions x in the nucleus probed by the deuteron. This trend is seen more 
clearly in Fig. [2] where JdA is plotted versus x A ™ 9 = ((pTi)e~^ + (pT2)e~^) / ^/saw 
for all pair selections in rj and pr- In the case of 2— >2 leading order (LO) processes, the 
variable x A ™ g is lower than x Au by the mean fragmentation fraction, (z), of the struck 
parton in the Au nucleus. Since x A ™ 9 is an entirely experimental defined quantity, it 
should be reproducable in any theoretical framework. 

3. Discussion 

In a leading order pQCD picture, the variable JdA is 

, = ° P d 7/°<iA „ fS(x a d ) ® f b Au (x Au ) ® d^ d g V(z c , z d ) 
dA (N coll ) a p p r/o- pp ~ (N coll ) ® f b M) ® ° ab ^ cd ® V(*e, *d) 

for partons a+b going to outgoing jets c+d, which then fragment to hadrons with 
longitudinal fractions z c , z d . In the above convolutions over p+p and d+A, most of the 
terms are expected to be roughly similar between p+p and d+Au except for the nuclear 
gluon pdf. Naively, JdA might be largely dominated by the modification to the nuclear 
gluon parton distribution function (pdf's), since most of the events with di-hadrons at 
forward rapidities consist of a high-x parton from the deuteron and a low-x gluon from 
the gold nucleus. Assuming this to be true then J dA ~ R^ u = Gau(x,Q 2 )/AG p (x,Q 2 ) 
In figure [2] the JdA values are overlaid with the EPS09 Rf u curves [15] . The Jja values 
for the peripheral bins are above the best fit EPS09, while the central bins are below. 
The EPS09 curves are taken largely from nuclear deep inelastic scattering and represent 
an averaged value of R^ u over all centralities. The Jja values for the most central bin 
at the lowest x are well below the EPS09 curves. This is quantitatively consistent with 
the expectations for the Color-Glass Condensate [3], which posits an extreme form of 
shadowing due to the onset of gluon saturation. 

If nature is kind and this data can be interpreted in terms of a simple LO 
pQCD picture, it may be possible to extract Rf u , which is extremely important for 
understanding the quark gluon plasma since it forms the baseline for production in 
heavy ion collisions. In addition, if the large suppression of Jja observed in central 
d+Au collisions is from gluon saturation, it may be possible to study the dependence 
of that saturation on the thickness of the nucleus. One possible test of whether these 
ideas are correct would be to use extractions of R^ u from this data to predict J/^f data 
in d+Au collisions from PHENIX. 
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